Evolution of a barotropic shear layer into elliptical vortices 
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Barotropic shear layers, on becoming unstable, produce the well known Kelvin-Helmholtz insta- 
bility (KH). The nonlinear manifestation of KH is usually in the form of spiral billows. However, 
a piecewise linear shear layer produces a different type of KH characterized by elliptical vortices 
of constant vorticity connected via thin braids. Using DNS and Contour Dynamics, we show that 
the interaction between two counter-propagating vorticity waves are solely responsible for this KH 
formation. We investigate the rotation and nutation of the elliptical vortex, and the oscillation and 
straining of the braids. Our analysis also provides possible explanation behind the formation and 
evolution of elliptical vortices appearing in geophysical and astrophysical flows, e.g. meddies, Jovian 
vortices, Neptune's Great Dark Spot, coherent vortices in the wind belts of Uranus, etc. 
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I. INTRODUCTION 

Barotropic shear layers are ubiquitous in the atmo- 
sphere and oceans. These layers can become hydrody- 
namically unstable, giving rise to an instability mecha- 
nism widely known as the Kelvin-Helmholtz instability 
(KH). The non-linear manifestation of KH is usually in 
the form of spiralling billows, breaking of which generates 
turbulence and subsequent mixing in geophysical flows. 

In theoretical and numerical studies, the hyperbolic 
tangent velocity profile (U(y) = tanh(y)) is often consid- 
ered to be the prototype of smooth barotropic shear lay- 
ers pQ . Initially interested in understanding the long time 
evolution of KH emanating from this profile, we have per- 
formed a moderate Reynolds number Direct Numerical 
Simulation (DNS); see Fig. [TJa). The flow re-laminarizes 
once the KH billow completely breaks down into small- 
est scales via turbulent processes. At this stage the shear 
layer is much sharper and around four times thicker than 
the initial tanh(y) profile. This four-fold increase in thick- 
ness also leads to an equal increase in the flow Reynolds 
number, making the flow more inviscid. Moreover, the 
shear layer profile at this stage closely matches the piece- 
wise linear profile (see Fig. [TJb,c)): 



U{y) = { 



i y> i 

y -i<y<l 

-i y<-i 



(l) 



Here U is the non-dimensional velocity of the shear flow. 
Eq. ([I]) now becomes the new base flow and serves as the 
initial condition for the subsequent instability processes. 

The linear stability analysis of the base flow given by 
Eq. ([!]) dates back to Lord Rayleigh [2]. The first non- 
linear analysis however was performed more than a cen- 
tury later. Using a boundary integral method known as 
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the Contour Dynamics, Pozrikidis and Higdon [3 showed 
that the piecewise linear shear profile evolves into nearly 
elliptical patches of constant vorticity - Kirchhoff vor- 
tices. We hypothesize that the initial shear layer pro- 
file determines the asymptotic form of the ensuing KH 
- smooth shear layers give rise to spiral billows, while 
sharp shear layers produce Kirchhoff vortices. 

The spiralling billow form of KH has been thoroughly 
investigated in the past. In fact, the spiral billow shape 
has become the signature of KH [3] ; but little is known 
about the non-linear evolution of the sharp shear layer. 
This is because the piecewise linear profile is usually con- 
sidered to be of little practical relevance. Its usage is 
therefore restricted to theoretical studies, mainly as a 
crude approximation of smooth shear layers [SJ [3] . How- 
ever, our DNS results in Fig. [TJa,b) indicate otherwise - 
the sharp profile can also occur in nature. The fact that 
this profile produces elliptical vortices similar to those ob- 
served in geophysical and astrophysical flows, e.g. med- 
dies in Atlantic ocean [7J, Neptune's Great Dark Spot 
[8], coherent vortices in the atmosphere of Jupiter and 
Uranus jS], has motivated us to investigate further. 



II. LINEAR THEORY 

In 1880, Lord Rayleigh [2] performed a linear stabil- 
ity analysis of the shear layer profile in Eq. ([!]) and 
showed it to be unstable for the range of wavenumbers 
< a < 0.64, the fastest growing mode being a cr i t = 0.4. 
In the conventional linear stability approach, infinitesi- 
mal wavelike perturbations are superimposed on a lami- 
nar background flow and an eigenvalue problem is solved 
to find the band of unstable wavenumbers 10J. This 
mathematical exercise, however, provides little insight 
into the underlying physical mechanism. In the past 50 
years there has been a continuous effort to provide a 
mechanistic picture of hydrodynamic instabilities, espe- 
cially homogeneous and stratified shear instabilities. An- 
alytical studies performed on the profile in Eq. ([I]) have 
shown that the interaction between two vorticity waves 



is responsible for the development of KH [5J ITlTfTo] . We 
refer to this wave interaction based interpretation of hy- 
drodynamic stability as the "Wave Interaction Theory" 
(WIT). 

KH can be understood in the light of WIT by referring 
to Eq. ([lj. The vorticity O = dU/dz is discontinuous 
at z = ±1, which allows each of these two locations to 
support a stable, progressive, interfacial wave called the 
vorticity wave (also known as the Rayleigh wave). In ro- 
tating frame, its analogue is the Rossby edge wave which 
exists at the discontinuities in Potential Vorticity. In gen- 
eral, a vorticity wave at the interface z = Zj has a phase 
speed c r (zj) given by 

c r (z j )=U(z j ) + c™ t (z j ) (2) 

The first component, U(zj), is the background velocity, 
while c*"*(zj) is the intrinsic phase speed defined as 

where [Ct] z . = Q(zj + e) — 0,(zj — e) is the jump in 
fl across Zj. For the profile in Eq. ([I]), the arrange- 
ment is such that c l r nt (zj) is always directed counter to 
U(zj). We define such waves as counter-propagating vor- 
ticity waves (CVWs) (analogous to counter-propagating 
Rossby waves in a rotating frame). The CVWs at z = ±1 
have [15] : 

c± = c r (±1) =±1T^ (4) 

The above equation implies that these two waves always 
travel in opposite directions. The wave at z = 1 is left 
moving while the one at z = — 1 moves to the right. 
Although each of them is stable, their interaction can 
lead to instability, producing KH. 
Consider a pair of CVWs: 

r] + = ao cos (ax — $) at z = 1 (5) 
rf~ = ao cos (ax) at z = — 1 (6) 

where the perturbation amplitude ao <C 1 and the phase 
shift $ G [— it, tt]. Modal growth occurs only when the 
two waves together behave like a normal mode. Two con- 
ditions need to be precisely satisfied for modal growth: 
(a) phase locking - the waves are stationary relative to 
each other, i.e. and c~ become zero after interac- 
tion, and (b) mutual growth - the phase shift between 
the waves is such that one makes the other grow ex- 
ponentially. Linear theory predicts this phase shift to 
be $ modal = cos- 1 {(1 - 2a)e 2a } [3. For a = a crU , 
the corresponding $™° t dQ ' = 0.353tt. Thus WIT provides 
the appropriate initial conditions (i.e. a cr i t and <& cr it) for 
modelling the non-linear evolution. 

Unfortunately, WIT itself is limited to the linear 
regime only. Current methodology does not allow a 
straightforward non-linear extension. This is a direct 
consequence of linearization which forces the vorticity 
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FIG. 1. (a) 3D DNS at Re = 2400 performed to capture the 
complete turbulent dissipation of a KH billow ensuing from 
an initial flow given by U(y) = tanh(y). False colour is added 
to aid visualization, (b) The dashed line represents the initial 
velocity profile, while the solid line shows the velocity profile 
once the flow re-laminarizes (t = t^). (c) The Magenta line is 
the continuous velocity profile obtained from Eq. (JSp , while 
the thick black line below it is the piecewise linear profile from 
Eq. |TJ. Both of them closely mimic the solid green line, (d) 
Linear vorticity waves (exaggerated) existing at the vorticity 
discontinuities. 

interfaces to become vortex sheets (H QTj [16]. There- 
fore it is questionable whether non-linear evolution KH 
can be understood in terms of two interacting non-linear 
vorticity waves. 

III. NON-LINEAR FORMULATION 

Over the last few decades, sophisticated computational 
techniques have been developed for precise understand- 
ing of turbulent processes. Two such techniques worth 
mentioning are Direct Numerical Simulation (DNS) and 
vortex methods. 

Many geophysical and astrophysical flows can be as- 
sumed incompressible and quasi-inviscid (involve high 
Reynolds number). In such flows, non- linear processes 
like chaos and turbulence are solely driven by vortex 
interactions [T7]. Vortex methods are especially useful 
under such circumstances; they numerically solve the in- 
viscid and incompressible Navier-Stokes equations (Euler 
equations). An example of one such 2D vortex method 
is Contour Dynamics [T5] , 

A. Contour Dynamics 

High Reynolds number flows have a tendency to de- 
velop Finite Area Vortex Regions (FAVR) with steep 
sides [IS]. If the flow is incompressible and 2D, the area 
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as well as the vorticity of an FAVR are conserved quan- 
tities. The contour (FAVR or "vortex patch" boundary) 
is basically a vorticity jump and a material surface [T7] . 
Substantial simplification is possible for constant vortic- 
ity patches; the governing 2D Euler equations can be 
reduced to a ID boundary integral. This provides a sig- 
nificant computational advantage of computing the evo- 
lution of an FAVR by only solving the contour motion. 
This methodology is known as the Contour Dynamics 
(CD). Numerical implementation of CD is done in a La- 
grangian framework by tracing the contour with a set of 
N marker points. 

Observing that a piecewise linear shear layer can be 
represented by a horizontally periodic patch of constant 
vorticity, Pozrikidis and Higdon [3] have used CD to sim- 
ulate the non-linear evolution process. The evolution of 
the i-th Lagrangian marker is given by [5] : 



G?Xj 

~dt 



47T 



In 



cosh ( aA 



cos 



dx (7) 



where fi = 1, x = [x, y] T , As- = x i - x , Ay- = y, - y 
and C is the contour around one patch. On perturbing 
the contour with sinusoidal disturbances, the shear layer 
rolls up producing nearly elliptical patches of constant 
vorticity [3] . Thus we infer that the contour is not just 
a material boundary, it is actually a pair of interacting 
CVWs. 

To elucidate this point further we refer to WIT which 
explains how the interaction between two CVWs, each 
existing at a vorticity discontinuity, solely governs the 
(linear) evolution of a piecewise linear shear layer. Like- 
wise, CD simulation has shown that the shear layer con- 
tour (like CVWs, the shear layer contour also exists at 
the vorticity discontinuities) is solely responsible for the 
(non-linear) evolution of a piecewise linear shear layer. 
It follows that the shear layer contour and the pair of 
CVWs are one and the same thing. This observation has 
another important implication - CD acts as a non-linear 
extension of WIT. Thus the shortcomings of WIT arising 
from linearization can be avoided, making it possible to 
understand the non-linear evolution of KH in terms of 
two interacting CVWs. 

We perform a CD simulation with initial condition 
based on WIT; we choose a — a crit and <I> = $™° t da/ , and 
initialize the problem using Eqs ([5])- (16]). This allows us to 
simulate the non-linear evolution of the fastest growing 
mode of KH arising in a piecewise linear shear layer. Con- 
sidering a domain one wavelength (\ cr u = 2Tr/a cr n = 
5ir) long, we solve Eq. ufh using central differencing for 
space derivatives and 4 order Runge-Kutta for time. 
Each wave is initially represented by 400 points. During 
its evolution process, an adaptive point insertion-deletion 
algorithm is used to check if the neighbouring points are 
within a desired distance. 

Pozrikidis and Higdon [3 studied the nonlinear evolu- 
tion of a piecewise linear shear layer subject to different 
types of initial disturbances, i.e. different combinations 



of wavenumbers (a) and phase shifts ($). For each ini- 
tial condition, the shear layer evolved into a correspond- 
ing nonlinear state. Since the wavenumber of maximum 
growth (a cr it) usually determines the non- linear state of 
a shear layer, it is unlikely to physically observe most of 
the non-linear states depicted in Pozrikidis and Higdon. 



B. Direct Numerical Simulation 

Previous studies have used CD mainly as a tool for 
qualitative understanding of problems involving inviscid 
vortical flows [19]. In order to demonstrate the quanti- 
tative capabilities of CD, we validate our CD simulation 
against a pseudo-spectral Direct Numerical Simulation 
(DNS). We also expect that DNS will provide a parallel 
understanding of the phenomena. The DNS code uses 
full Fourier transform in the horizontal direction, and 
half-range sine or cosine Fourier transform in the vertical 
direction in order to convert the set of partial differential 
equations (Navier-Stokes equations) into ordinary differ- 
ential equations. Time integration is performed using a 
third-order Adams-Bashforth method. Detailed descrip- 
tion of this code can be found in the work of Winters et 
al. [20]. 

We consider a domain of length X cr it in the horizon- 
tal direction and nine times the initial shear layer thick- 
ness in the vertical direction. The horizontal boundary 
condition is periodic while the vertical boundary condi- 
tion is no-flux free-slip. We perform a 2D simulation at 
Reynolds number Re = 10, 000 (Re — 1/z/, where v is the 
fluid viscosity). At such high i?e, DNS provides a good 
representation of CD. Our simulation is resolved using 
2880 x 3456 points, which is fine enough to simulate the 
smallest scales of the 2D flow. Since non-differentiable 
profiles like Eq. cause numerical issues (Gibbs phe- 
nomena), we use a smooth velocity profile that resembles 
the piecewise linear profile accurately; see Fig.[TJc). This 
velocity profile is derived from a vorticity distribution 
having the form 



n( y ) = 



1 



1 — tanh 



y 2 -i 



(8) 



Integrating Eq. (|8]), the velocity profile is obtained di- 
rectly: U = j fldy. The linear stability characteristic of 
this profile matches almost exactly with the piecewise 
profile, and has the same a cr n. Equating the total cir- 
culation of DNS with CD yields e = 0.100. The vorticity 
field in DNS is perturbed to match the initial wave am- 
plitude growth in CD. 



IV. RESULTS AND DISCUSSION 

A. Pre-saturation and Saturation Phases 

The non-linear evolution of KH is illustrated in Fig. [2] 
It shows the vorticity field from DNS and contour lines 
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FIG. 2. Time evolution of Kelvin- Helmholtz instability - comparison between DNS (top) and CD (bottom). 




shear layer between the crest of the lower wave and the 
trough of the upper wave (see Fig. [ljd)) gives rise to the 
elliptical core, the initial length of which is given by 



FIG. 3. Temporal variation of the wave amplitude a. The 
blue line is the prediction from linear theory, the black line 
corresponds to CD while the magenta line represents DNS re- 
sults. The inset shows the variation of the ellipse aspect ratio 
r (green line) and the angular rotation rate uj (red line) with 
time, obtained by solving Eqs. (JToJl - (jXTJ . The black markers 
indicate the data points measured from DNS. 



from CD. Vorticity of the region enclosed by the con- 
tour lines is conserved in CD. However the presence of 
viscosity makes conservation of vorticity invalid in DNS. 
The implementation of high Re DNS minimizes the vis- 
cous effects, making DNS comparable to CD. The basic 
premise behind our simulations is the conservation of to- 
tal circulation T — HA (where the vorticity ft — 1 and A 
is the shear layer area), which comes from Kelvin's Cir- 
culation theorem |T7] . Its corollary is the conservation 
of shear layer area - a quantity that remains fixed at its 
initial value A = 2X cr i t . 

Fig. [3] shows the time evolution of the wave amplitude, 
a. The maximum shear layer thickness, H — 2(1 + a) 
evolves in a fashion similar to the wave amplitude a. We 
find the growth to be exponential, at least for t ;$ 20. 
CVW interaction causes the shear layer to grow non- 
linearly. This phenomenon leads to the roll-up and for- 
mation of the elliptical core vortex (Kirchhoff vortex). 
The evolution process is shown in Fig. [2] The part of the 



fomodal 

crit 



(9) 



The flow saturates (i.e. the amplitude reaches a max- 
ima) at t sat — 34. For t > t sat , approximately 80% of 
F is concentrated in the core. H also reaches a maxima 
at saturation, and has the value H max — 8.7. The fully 
formed elliptical cores are connected by thin filaments of 
fluid called braids. These braids wrap around the rotat- 
ing cores, producing a complex spiral like structure, as 
shown in Fig. [5] 



B. Early Post-saturation Phase 

After saturation, the core rotates with an angular ve- 
locity u>, causing the wave amplitude, a, to oscillate with 
a time period T amp ss 13; see Fig. [3] The core also nu- 
tates, i.e. the core aspect ratio r (defined as the ratio 
between the ellipse major axis and the minor axis) un- 
dergoes a periodic oscillation. This phenomenon is called 
nutation and is apparent from both Figs. [2] and [3] 



1. Nutation 

To better understand the nutation process, we consider 
the simple model proposed by Kida [3T]. An isolated 
Kirchhoff vortex rotates in the presence of a constant 
background strain-rate 7. The velocity field associated 
with this strain-rate is given by u s — ja, w s = —7^ where 
a and £ are the principal axes with the origin at the cen- 
tre of the ellipse. In our case, this velocity field mimics 
the leading order straining effect induced by the rota- 
tion of other Kirchhoff vortices. Note that the periodic 
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boundary condition takes into account the effects of the 
other Kirchhoff vortices. Let the clockwise angle between 
a and the ellipse major axis be 9 at any instant. Then 9 
and r evolve as follows [2"T] : 



UJ 



(19 



"7 



1 



sin (26») 



Sir 



(r + iy 



dr 
dt 



2^r cos (29) 



(10) 
(11) 



Eq. (Ill implies that the nutation is caused by strain. It 
also reveals that r reaches maxima at 9 = ±7r/4. Si- 
multaneously, Fig. [2] shows that the core nutates with 
a maximum value of r along the x axis and a minimum 
along the y axis. Therefore the a axis must make an angle 
7r/4 with the x axis. The angle made by the braid with 
the x axis at the stagnation point(s) is also 7r/4. This is 
because the braid aligns itself with the streamlines. 

To answer why the a axis makes an angle of 7r/4 with 
the x axis, we consider an ideal problem where an infinite 
number of Kirchhoff vortices, each of circulation T core = 
T (note T = 2X crit ), are placed along the x axis with a 
constant spacing X cr it between their centres. The effect 
of a vortex patch at a point outside its own boundary is 
exactly the same as that of a point vortex of equivalent 
strength placed at the center of the patch [5] . Therefore 
we replace all the Kirchhoff vortices with point vortices of 
strength r core . This provides a simplistic understanding 
of the mechanism by which the rotation of distant vortex 
patches strain a given patch. We find this ideal strain- 
rate to be 



E < 

— og , n ^ 



0.067 



(12) 



The principal axes of the strain field produced by this 
infinite array of point vortices make angles of ±7r/4 with 
the x axis. Hence, this ideal strain field and the strain 
field of our actual problem have the same orientation. 
Before comparing the magnitudes of these two fields, 
it is important to note that the presence of the braids 
complicate the actual problem by making the strain-rate 
magnitude vary spatially. We simplify the analysis by 
assuming a strain field of constant magnitude acting on 
the elliptical core, thereby reducing the problem to the 
Kida problem described by Eqs. ([H^-([TT]). DNS is used 



to supplement the analysis by providing the values of r 
and 9 wherever necessary. By applying this methodol- 
ogy, the magnitude of the actual strain-rate is found to 
be 7 = 0.073, which is quite close to the ideal value of 
0.067 obtained from Eq. Q . 

We also capture the evolution of r and uj by solving 
Eqs. ( X0 )-( 11 1; see Fig. [3] The initial values are obtained 
from DNS, and 7 = 0.073. The spike in uj at t « 45 
is caused by r — > 1. The data points measured from 
DNS when compared with the Kida theory reveal that 
the latter can provide encouraging predictions. The nu- 
tation period is found to be T nut w 13, while the period 
of core rotation is T core — 2tt/uj w 26 (overbar denotes 



average). T core ss 2T nut is because one full rotation cor- 
responds to passing the coordinate axes twice. Likewise, 
T core ss 2T amp because the braids are connected to the 
two ends of the core. 



2. Small Length Scale Production 

The smallest length scales are found to occur in the 
braid region adjacent to the core; see Fig. [2] This is due 
to the straining effect of the rotating elliptical core which 
causes the braid region in its vicinity to thin exponen- 
tially fast. In real flows, when a fluid element becomes 
sufficiently thin, the balance between the strain-rate and 
the viscous dissipation determines the small length scales. 
The order of magnitude of the core rotation induced 
strain-rate, jiocaU can b e obtained by replacing the core 
with a point vortex of equivalent strength and located at 
the ellipse centre: 



'Ylocal 



r 1 

1 core ± 



2irl 2 



(13) 



where I is the characteristic length of the core. No- 
tice that this local strain-rate is one order of magnitude 
greater than the background strain-rate 7 or 7'. 

The smallest length scale appearing in a 2D turbulent 
flow is L2D ~ -Re -1 / 2 [22]. This length scale is a 2D ana- 
logue of the Taylor microscale occurring in 3D turbulent 
flows. In order to estimate the time when Lid appears in 
our flow, we formulate a braid evolution equation similar 
to Eq. (2.8) of Corcos and Sherman 23 : 



S 2 (t*)=S 2 (0)e 



-2 7loc t* 



2j ioc Re 



1 



,-2ji oc t' 



(14) 



where S (£*) is the braid thickness adjacent to the core 
at time t* = t — t sat . We estimate 6 (0) from DNS and 



solve Eq. ( 14 ) . We find that L2D appears soon after the 
saturation, around t* ~ 4, implying that 2D transitional 
flows like KH can give rise to "turbulent" features at a 
very early stage. 



C. Late Post-saturation Phase 

For t 50, the core surface develops progressive vor- 
ticity waves. These are called the Kelvin m- waves or 
the Love to- waves [24] . where to is the eigenmode. In 
the absence of strain, r < 3 is the condition for stabil- 
ity of a Kirchhoff vortex [T7]. Although r satisfies this 
condition in our case (refer to the inset in Fig. [3|, the 
presence of strain adds instability [25]. Finite ampli- 
tude to = 4 waves give rise to winding filaments, see Fig. 
|4j This is a feature of Kirchhoff vortices in background 
shear or strain, which has been thoroughly investigated 
by Dritschel [25]. Background shear or strain is how- 
ever not necessary to produce winding filaments. Similar 
winding features are also seen in large aspect ratio Kirch- 
hoff vortices even in the absence of shear or strain 1241. 
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t = 67 




FIG. 4. Formation of winding filaments around the elliptical 
vortex during the late post-saturation phase. 



V. PRACTICAL IMPLICATIONS 

In this Section we discuss several practical aspects of 
our work. The elliptical core shown in our simulations is a 
common feature in geophysical and astrophysical flows, 
e.g. coherent vortices in the wind belts of Uranus [3], 
Mediterranean eddies or meddies [7J, vortices in Jupiter 
[5J , the Great Dark Spot of Neptune (which has dis- 
appeared now) [8], etc. 

The wind belts of Uranus often produce an array of co- 
herent elliptical vortices, see Fig. 3 of Liu and Schneider 
[9J . The probable mechanism behind the formation of the 
Uranian vortices, as suggested by our CD and DNS simu- 
lations, is the interaction of counter-propagating Rossby 
waves. Two dimensional CD and DNS modelling with 
streamwise periodic boundary condition can be justified 
in this context which implies that our simulations can 
provide a good description of the Uranian vortex dynam- 
ics. However, an added background shear present in the 
wind belts of Uranus slightly complicates the problem. 
The effect of this added shear can be incorporated in 
CD, DNS and the Kida model, making the subsequent 
analysis a straightforward extension of our present study. 

The Great Red Spot and the White Oval BC of Jupiter, 
and the Great Dark Spot of Neptune can be well ap- 
proximated as Kirchhoff vortices in a uniform back- 
ground shear [8]. The rotation and nutation of these 
vortices were analyzed by Polvani et al. [5J using the 
Kida model. We hypothesize that Jovian and Neptu- 
nian vortices are produced by the interaction of counter- 
propagating Rossby waves in a fashion which is similar 
to the Uranian vortices. 

Elliptical vortices known as the meddies are found 
in the region where Mediterranean sea meets the At- 
lantic ocean. The genesis of meddies is highly speculative; 
one possible explanation is the interaction of counter- 



propagating Rossby waves. This is corroborated by the 
observations near Portimao Canyon and Cape St. Vin- 
cent [7J. It is possible that the meddies, like Uranian 
and other giant planetary vortices, are also a non-linear 
manifestation of KH ensuing from sharp shear layers. 

The average diameter of a meddy is 25 — 100 km which 
is approximately half of the most unstable wavelength of 
the Rossby edge waves [7J . A simple visual inspection of 
Fig. [^corroborates this observation. During its evolution 
process, meddies are found to produce winding filaments 
[2"5] , similar to what we observe in Fig. |4j However for 
proper understanding of meddy evolution, density strat- 
ification, three dimensionality and aperiodic streamwise 
boundary condition should be taken into consideration. 

VI. CONCLUSION 

A piecewise linear shear layer becomes unstable and 
evolves into a series of elliptical vortices of constant 
vorticity (Kirchhoff vortices) connected by thin braids. 
The interaction between two counter-propagating vortic- 
ity waves is the driving mechanism behind this instability 
process. Although this fact was known previously, lin- 
earized approximations forced the analysis to be valid 
only in the linear regime. By finding and exploiting 
the link between two quite different theories, namely the 
Wave Interaction Theory and the Contour Dynamics, we 
are able to extend the analysis to the fully non-linear 
domain. 

The production of Kirchhoff vortices clearly indicates 
that KH arising from a piecewise linear shear layer is very 
different from the classical spiralling billow type KH en- 
suing from a smooth shear layer. The characteristics of 
this little known KH have been investigated. The ro- 
tation and nutation of the Kirchhoff vortices are found 
to be consistent with the predictions of Kida. The time 
period of rotation of the vortices is found to be twice 
the period of nutation and the period of maximum shear 
layer height oscillation. The braids connecting the Kirch- 
hoff vortices rapidly thin to a length scale which is the 
2D equivalent of Taylor microscale. 

Elliptical vortical structures, similar to those found in 
our simulations, are quite common in nature, especially 
in regions with sharp shear. Examples of such vortices 
include meddies, Jovian vortices, Neptune's Great Dark 
Spot, coherent vortices in the wind belts of Uranus, etc. 
Our analysis can help understanding their formation and 
evolution process, motivating further investigation. 
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